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Abstract 

We make a theoretical study of the s-wave interaction of the nonet of vector 
mesons with the octet of pseudoscalar mesons starting from a chiral invariant La- 
grangian and implementing unitarity in coupled channels. By looking for poles in 
the unphysical Riemann sheets of the unitarized scattering amplitudes, we get two 
octets and one singlet of axial-vector dynamically generated resonances. The poles 
found can be associated to most of the low lying axial-vector resonances quoted in 
the Particle Data Book; 6i(1235), /ii(1170), /ii(1380), ai(1260), /i(1285) and two 
poles to the Rri(1270) resonance. We evaluate the couplings of the resonances to the 
VP states and the partial decay widths in order to reinforce the arguments in the 
discussion. 


1 Introduction 

The realization that QCD at low energies can be studied by means of effective chiral 
Lagrangians using helds associated to the observables mesons and baryons [1-4] has brought 
a substantial progress in hadronic physics by means of chiral perturbation theory. A further 
useful step in this direction has been given with the introduction of unitary techniques 
which have allowed to extend to higher energies the predictions of chiral perturbation 
theory, as well as to tackle new problems barred to a pure perturbative expansion [5- 
10]. The unitary extensions of chiral perturbation theory, t/yPT, have brought new light 
in the study of the meson meson and meson baryon interaction and have shown that 
some well known resonances qualify as dynamically generated, or in simpler words, they 
are quasibound states of meson meson or meson baryon. This is the case of the low 
lying scalar mesons a, /o(980), ao(980), fi:(900) [8-12], which appear from the interaction 
of pseudoscalar mesons. Another case of successful application of these chiral unitary 
techniques is the interaction of mesons with baryons [6,7,13-16] showing that the A(1405) 
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and the A^*(1535) were dynamically generated resonances. A more systematic stndy of 
these latter interaction has shown that there are two octets and one singlet of resonances 
from the interaction of the octet of psendoscalar mesons with the octet of stable baryons 
[17,18]. Work along these lines has continned by stndying the interaction of the octet 
of psendoscalar mesons with the decnplet of baryons [19, 20] which has also led to the 
generation of many known resonances, like the A^*(1520) and the A(1520). 

The stndies of the meson meson and meson baryon interaction along these lines have 
also shown that some mesons or baryons are not dynamically generated, they are not 
conseqnence of the interaction between the meson or baryon components and they qnalify 
better as gennine, or preexistent states, a word that can be snbstantiated as that they 
wonld remain in the limit of large where the loops of intermediate states vanish. This 
is the case of the vector mesons of the p octet [10,21]. 

A snggestive step forward in the direction of stndying the meson meson interaction 
along these lines would be to study the interaction of the octet of psendoscalar mesons 
with the octet of vector mesons (nonet including 0 and uj with the standard mixing). 
Work in this direction has already been done in [22] and, interestingly, most of the low 
lying axial-vector mesons appear as dynamically generated resonances. If this were the 
case, this would have many practical as well as conceptual repercussions and, hence, extra 
efforts to corroborate these hndings, looking also for uncertainties, are called for. In [22] 
the search for resonances was done by looking at the speed plots of the physical amplitudes. 
A search for poles of the amplitudes in the unphysical Riemann sheets in the complex plane 
is a more powerful tool to investigate resonances an a thorough study along these lines is 
called for. The present work goes in this last direction and we have done a thorough work 
investigating the following points: 1) The poles of the amplitudes have been systematically 
searched for in the complex energy plane following the trajectories in terms of an SU(3) 
breaking parameter. This has allowed us to make an SU(3) study of the problem as well as 
the effects of SU(3) symmetry breaking, very useful to understand the meaning of the poles 
and their number. 2) The study of the amplitudes around the poles has also allowed us to 
determine the coupling of the resonances to the different coupled channels and from there 
the partial decay widths into the different channels. This study has been useful to make 
an association of the resonances found with those of the Particle Data Group (PDG) [23], 
(see Table Q). 3) We have changed the association of the poles to the known resonances in 
some cases with respect to [22], in particular in the case of the iPi(1270) resonance, which 
we claim comes from two distinct poles with very different properties. This has some 
practical consequences in the partial decay widths and the dominance of certain channels 
in different reactions, which would help clarifying the puzzle of these resonances. 4) At 
the same time, and in order to help new researchers in the held, we have substantially 
simplihed the formalism of [22], making it more amenable and transparent. 

The results obtained here support the bulk of the claims made in [22], but the larger 
amount of information on these resonance obtained here brings new elements that set 
on hrmer grounds the association of these resonances with the experimental ones. This 
substantiates the idea that the low lying axial-vector resonances are dynamically generated, 
with the exception of the higher mass ones, the /i(1420) and the iPi(1400), which do not 


2 


fit easily in our scheme, while at the same time we suggest that the ii'i(1270) corresponds 
actually to two poles, which would have many experimental repercussions. 


JPC 

/ = 1 

1 = 0 

/ = 1/2 

1+- 

6i(1235) 

/ti(1170), /ti(1380) 

Kib 

Ty /i Ty (^ /inn\ 

1++ 

ai(1260) 

/i(1285), /i(1420) 

itifl4UUj 

P-iA 


Table 1: I"*" resonances in the Particle Data Book [23]. The ii'i(1270) and 77i(1400) are 
assumed to be a mixture of Kia and Kib- 


2 Pseudoscalar-vector meson interaction 

2.1 Tree level potential 

There is not a unique formulation to treat the vector mesons in an effective theory. The 
ambiguity comes from the freedom in considering the nature of the vector meson to be or 
not a gauge boson of a certain symmetry, and from the election of their transformation 
properties under a certain realization of chiral symmetry. For a review see, for example. 
Ref. [24]. Despite the differences between the treatments of vector mesons, the equivalence 
between the different approaches can be shown at lowest order [24]. Considering the 
vector mesons as helds transforming homogeneously under the nonlinear realization of 
chiral symmetry, the interaction of two vector and two pseudoscalar mesons at lowest order 
in the pseudoscalar helds can be obtained from the following interaction Lagrangian [24]^ 

Cl = —Tr {(V^K - - V"R^)} , (1) 

where Tr means SU{3) trace and is the covariant derivative dehned as 

V^K = 9^K + [r^,K], (2) 

where [, ] stands for commutator and T^ is the vector current 

CfM = + udf,u^) (3) 

with 

u^ = U = (4) 

In the previous equations / = 92 MeV is the pion decay constant and P and V are 
the SU{3) matrices containing the octet of pseudoscalar and the nonet of vector mesons 
respectively: 

^Note the different factor —1/4 instead of —1/2 in [24] to agree with our normalization of the fields. 
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The Lagrangian of Eq. ([Q) is invariant under the chiral transformations SU{3)l ® 
SU{3)fi, since transforms as [25] 

^ ( 6 ) 

We are interested in the two-vector-two-pseudoscalar amplitudes. Hence, expanding the 
Lagrangian of Eq. (HD up to two pseudoscalar meson fields we find 

Cvvpp = -^TrilV^, cfV„]\P,d^P ]), ( 7 ) 

which would account for the Weinberg-Tomozawa interaction for the VP ^ VP process 
[22,24,26,27]. 

Note that in Eq. © in the pseudoscalar octet we are only considering the rjg = rj and 
not the 7]'. The inclusion of the rj' effects in strong interactions can be accommodated 
in the higher order Lagrangians [28]. Since the meson decay constant, /, is different for 
different mesons, one could also use different values of / for the different pseudoscalars, as 
done in [13,16]. We shall comment on the results obtained when we take this into account. 
In the vector meson multiplet we have assumed ideal uji — cug mixing: 

(f) = cui/a/S — ctJg2/3, u! = u)i\J‘2/3 -\- a;g/\/3. (8) 

(Throughout the work we will use the phase convention | 7 r+ >= —11 -|- 1 >, Ip"*" >= 
— 11 -I- 1 >, \K~ >= —|l/2 — 1/2 > and \K*~ >= —|l/2 — 1/2 > corresponding I/J 3 > 
isospin states). 

From the Lagrangian of Eq. © one obtains the full amplitude which we deduce in 
Appendix I, where we also make the projection over s-wave, which leads to 


r«(s)= 


6 • 6 

r 

8/2 


3s - (M" + + M'" + m'^) - -(M" - - m'") 

s 


(9) 


where e(e') stands for the polarization four-vector of the incoming(outgoing) vector meson. 
The masses M(M'), m{m') correspond to the initial(final) vector mesons and initial(final) 
pseudoscalar mesons respectively, and we use an averaged value for each isospin multiplet. 
The indices i and j represent the initial and final VP states respectively. 

We can identify the VP channels by its strangeness (S') and isospin (J), (S', I) = (1,1/2), 
(0,0) and (0,1). There are other possible (S',/) combinations but since, advancing some 
results, we will not find poles there, we will not consider them in the discussion. For the 
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(0, 0) channel the allowed isospin channels are K*K, (prj, ujrj, pn and K*K, for the (0,1) 
channels they are K*K, 07r, utt, pp, pvr and K*K, and for (1,1/2) we have (pK, uK, 
pK, K*p and K*n. Note that for the (0, 0) and (0,1) cases, the isospin states have well- 
dehned G—parity^ except the K*K and K*K states, but the combinations 1/\/2{\K*K > 
±\K*K >) are actually G—parity eigenstates with eigenvalues ±. 

In Tables 121 El and m we show the Cij coefficients in isospin base for (5,1) = (1,1/2), 
(0, 0) and (0,1), indicating also the G—parity in the (0, 0) and (0,1) cases. 
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Table 2: Cij coefficients in isospin base for S' = 1, J = |. 
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Table 3: Cij coefficients in isospin base for S' = 0, J = 0. The hrst column indicates the 
G—parity. 


Let us now discuss an interesting consequence in the sign and strength of the potential 
obtained, which can give us an indication about whether this Lagrangian can produce 
pseudo-bound states with a suitable unitarization procedure. The interaction of two octets 
gives the following decomposition in irreducible representations of SU{3): 

8 ® 8 = 1 © 8a © 8, © 10 © 10 © 27. (10) 

Note that although we shall work with a nonet of vector mesons (including the p and ui) 
the singlet component does not lead to an interaction term from Eq. O- 

^Recall that the G—parity operation can be defined through its action on an eigenstate of hypercharge 
(y), isospin (J), and third isospin projection {I 3 ) as G\Y,I,l 3 >= ~ >> with 77 being 

the charge conjugation of a neutral non-strange member of the SU(S) family. 
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Table 4: Cij coefficients in isospin base for S' = 0, / = 1. The first colnmn indicates the 
G—parity. 


The coefficients Cij of Eq. m can be expressed in any desired base: charge, isospin, 
SU(3), etc. Since the Sf/(3) matrices are given in terms of the physical charge fields, the 
most easy way to express the Cij coefficients is in this base, but with the use of the SU (3) 
Clebsch-Gordan coefficients it is straightforward to obtain the coefficients C in the SU (3) 
base, a. They give 

Cafs = diag{—6,—3,—3,0,0,2} (11) 

in the order of 1, 8a, 8^, 10, 10 and 27. As we shall see, a minus sign in a coefficient of 
Eq. (fTTjl implies an attractive potential, which is needed to have a bound state. Therefore, 
in the SU{3) limit, we should expect attraction in the singlet and the two octets, no 
interaction in the decuplets and repulsion in the 27—plet. Therefore, a priori, one could 
expect, after the unitarization procedure that we will explain below, two octets and one 
singlet of dynamically generated axial-vector (J^ = I’*') resonances. In addition, in the 
exact SU{3) symmetric case the two octets would be degenerate. 

2.2 Unitarization procedure 

In the literature several unitarization procedures have been used to obtain a scattering 
matrix fulhlling exact unitarity in coupled channels, like the Inverse Amplitude Method 
[5,9] or the N/D method [10]. In this latter work the equivalence with the Bethe-Salpeter 
equation used in [8] was established. 

In the present work we make use of the Bethe-Salpeter equation to resum the diagram¬ 
matic series expressed in Fig. [T] 

X '--x + '■( + 

Figure 1: Diagrammatic representation of the resummation of loops in the unitarization 
procedure. 
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There are some subtle differences with respect, for instance, to the pseudoscalar- 
pseudoscalar case, coming from the polarization vectors appearing in the potential and 
the particular form of the vector meson propagator in the loop. For the sake of clarity 
in the exposition we have relegated the explanation to Appendix II. From the reasons ex¬ 
plained in Appendix II we can do the evaluation of the scattering matrix for transverse 
polarization modes of the external vector mesons, which leads to the following unitarized 
amplitude: 


T =[l + VG]-\-V)e-e’, (12) 

where G = G(1 -|- 3 ]^) is a diagonal matrix with the /—th element, Gi, being the two 
meson loop function containing a vector and a pseudoscalar meson: 


GiiVs) = i j 


dS 


(27r)^ (P — qY — Mf + ie — rrif + ie'' 

with P the total incident momentum, which in the center of mass frame is (a/s, 0, 0, 0). 

The structure of Eq. oa explains why a minus sign in the Gaq coefficients of Eq. CD 
together with Eq. implies attraction, (ReG is negative in the region of relevance and 
1 + VG can lead to poles with V positive). 

In the dimensional regularization scheme the loop function of Eq. m gives 


Gii^Ts) = 




a(/i) -|- In 


Ml 

/i2 


mf — M? -|- s , mf 

-\ - - - - -In —- 

^ 2s Mf 


(14) 


+ 


qi 


s 


ln(s - {Mf - mf) + 2qi^/s) + ln(s -F {Mf - mf) + 2qiy/s) 


— ln(s — {Mf — mf) — 2qiy/s) — ln(s -|- {Mf — mf) — 2qiy/s) — 27ri 


where /r is the scale of dimensional regularization. Changes in the scale are reabsorbed in 
the subtraction constant a(/i), so that the results remain scale independent. In Eq. (HH), 
qi denotes the three-momentum of the vector or pseudoscalar meson in the center of mass 
frame and is given by 

qi = - {Mi + mi)"^]^ - {Mi - m/)^] (15) 

where Mi and mi are the masses of the vector and pseudoscalar mesons respectively. For the 
evaluation of the loop in the physical, or Erst, Riemann sheet one has to take the solution 
for the square root with Im(g/ positive. Note that in Eq. (ITT|l there is an ambiguity in the 
imaginary part of the In function coming from their multivaluedness. This ambiguity can 
be removed, for a generally complex ^/s, by comparing the result with the one obtained 
numerically by regularizing Eq. m by means of a cutoff of a natural size, of the order 
of 1000 MeV. By doing this, we have checked that the prescription for the In which gives 
a result in accordance with the cutoff method is to use the In with its argument defined 
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with a phase from —vr to vr with a cut in the negative real axis. On the other hand, this 
comparison with the cutoff method allows us also to determine the subtraction constant 
which turn out to be of ~ —1.85. 

In [22] a different regularization procedure is used by choosing the G function to vanish 
for ^/s equal to the mass of the vector meson. A similar choice, making the G function 
vanish at the mass of the baryon is shown to lead to realistic results in the meson baryon 
interaction case [18]. In the present work we will use both approaches in order to have an 
idea of the theoretical uncertainties in the results. 

On the other hand, note that we do not include the width of the vector mesons in their 
propagators in Eq. (USD. We shall take this into account in a coming section. 

2.3 Comparison to perturbative expansion 

This unitarization procedure can be understood as an analytical extrapolation of perturba¬ 
tion theory to higher energies in the same way as is done for ordinary chiral perturbation 
theory with pseudoscalar meson-meson or meson-baryon interaction [8,11,13,15]. Perturba¬ 
tion theory with the Lagrangian of Eq. dZD would proceed like ordinary chiral perturbation 
theory for pseudoscalar meson-meson and meson-baryon interaction, with loop divergences 
canceled with higher order Lagrangians (see [24]) which are ordered in terms of increasing 
number of derivatives in the helds. The expansion appears as a power series of the mo¬ 
mentum over one scale (expansion parameter), ~ dvr/ ~ 1.2 GeV (in the chiral limit 
of pseudoscalar masses going to zero). At one loop level one would need the next order 
Lagrangian to reabsorb the loop divergences and one would have the direct s-channel loop 
as well as the crossed loop term. The unitary amplitude should match the perturbative 
expansion at low energies, but our procedure only provides the s-channel loop and fur¬ 
thermore does not use information of higher order Lagrangians. The philosophy behind 
this is that the contribution of crossed loops generates a smooth energy dependence in the 
energy region of our concern and can be reabsorb in the subtraction constant, a, of the 
loop function G [10]. Similarly, one is also assuming that the one loop calculation with 
the lowest order Lagrangian, with the use of a G function with natural size subtraction 
constant, can account for the effect of the higher order Lagrangians. This is what char¬ 
acterizes the dynamically generated resonances, in contrast to cases like the p meson (a 
genuine resonance of basically two constituent quarks) which require the explicit use of a 
higher order Lagrangian [10]. 

In order to illustrate the previous discussion, we shall now compare results for the ampli¬ 
tudes obtained with the unitarity amplitude of Eq. m and its perturbative expansion up 
to two loops. We choose as an example the K*7i — K*tt in (5, 1) = (1,1/2) and ojtt —»• cutt, 
(S', I) = (0,1). First we show, in Fig. El results for the K*7i —> K*7i amplitude using only 
one channel and in the chiral limit (m^r = 0). In Fig. El we can see the modulus 

squared of the amplitudes calculated with the approximations (—E), {—V) + (—E)G(—E), 
{—V) + {—V)G{—V) + {—V)G{—V)G{—V) and the unitary amplitude. We can see that 
for low momenta there is a nice convergence of the perturbative series to the unitary result 
up to about 150 MeV/c. However, the unitary amplitude has a resonant structure with a 



Figure 2: Comparison of the perturbative expansion to the unitarized amplitude for K*^ 
K*tt with only one channel and in the chiral limit (m^ = 0). All the lines represent 

the modulus square. Dotted line: tree level potential (—D); dashed-dotted line: expansion 
of T at one loop, {—V) -f {—V)G{—V)] dashed line: expansion of T at two loops, {—V) + 
{—V)G{—V) + {—V)G{—V)G{—V)-, solid line: unitarized amplitude, T. 


peak around 250 MeV/c where the perturbative expansion is seen to fail drastically. Note 
that the lowest order {—V) (in the chiral limit) is of order q, as can be seen from Eq. 
or expanding Eq. ©• 

The real case has hnite pseudoscalar masses and coupled channels ujK, pK, K*ri, 
K*tt) and one has different thresholds for different channels. This case is shown in Fig. |21 
This fact changes the behavior of the amplitudes close to the K*^ threshold, where they no 
longer vanish, and, although the perturbative expansion is seen to converge to the unitary 
amplitude, the convergence is now slower. Once more we see a lack of convergence when 
we get close to the resonance peak around 170 MeV/c. 

In Fig. m we show the same amplitudes as before for the ujtt —>■ cutt cases. Now, the tree 
level amplitude is zero, see Table IH and the hrst resonance appears around 350 MeV/c. 
We see that the convergence of the perturbative expansion extends up to higher momenta 
than in the previous case, up to about 250 MeV/c. A combined conclusion from these 
two cases is that the perturbative expansion works well as far as we are not close to 
the lowest resonance, and of course, as far as q is small compared to A^. This reflects the 
mathematical theorem of complex variable that a series converges up to the hrst singularity, 
in our case the poles in the complex plane associated to the resonances. It is worth noting 
that what limits the convergence of the series expansion is the appearance of the lowest 
resonance, which in the present work appears at momenta much smaller than A^. The 
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Figure 3: Same as Fig.|21but with all the allowed VP channels and with finite pseudoscalar 
masses 


unitary amplitude has no problems of convergence, based as it is in the N/D method, 
as far as one used a proper interaction and included all relevant coupled channels. Up 
to energies of \/s ~ 1.5 GeV we consider that we are taking into account the relevant 
channels. This is about a momentum of about 500 MeV/c for th pvr channel, which has 
the highest momentum. This momentum is still smaller than the 600 MeV/c of the tttt 
interaction, where the chiral unitary approach was shown to be rather successful in [10]. 
Nevertheless we shall comment in the results sections on possible consequences from the 
inclusion of channels other than VP. 

Contrary to the perturbative expansion, the unitary amplitude allows for poles corre¬ 
sponding to bound states or resonances. In the next section we address this issue. 


3 Search for poles 

3.1 Unphysical Riemann sheets 

The association of physical resonances to poles of the scattering matrix in unphysical 
Riemann sheets is a very powerful tool to identify the resonances. The results of the 
scattering theory say that bound states reflect as a pole for Im(g) > 0 and Re(g) = 0, with 
q the momentum qi of Eq. (USD, i.e., in the real s axis below the lowest threshold. The 
resonances can appear only for Im(g) < 0 which means s with an argument larger than 
2n and Re(s) above the lowest threshold. This is what we will call second Riemann sheet 
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Figure 4: Same as Fig. IHlfor the ujn —^ (uvr case 


(R2) for the function T for the variable s. If these poles are not very far from the real 
axis they occur in -y/i = (M ± iF/2) = with M and F the mass and width of the 
resonance respectively. Of course the only meaningful physical quantity is the value of the 
amplitudes for real s/s, i.e., the reflexion of the pole on the real axis. Therefore, only poles 
not very far from the real axis would be easily identihed experimentally as a resonance. 

The effect of passing s to R2 has consequences only for the G functions. To evaluate 
G in R2 we can use the Schwartz reflexion theorem which states that if a function f{z) is 
analytic in a region of the complex plane including a portion of the real axis in which / 
is real, then [f{z*)]* = f{z). The loop function Gi satishes these conditions, therefore, for 
Re(yi) > mi + Ml we have 


G/(a/s - ie) = [G/(a/s + ie)]* = Gi{y/s + ie) - i2lmGi{y/s + ie). 
Since the beginning of R2 is equal to the end of R1 we have 

Gj\\/s + ie) = Gj (a/s — ie) = G\ (a/s + ie) — i2\mG{ (a/s + ie). 


(16) 


(17) 


where the superindices I and II refer to R1 and R2 respectively. 

The imaginary part of the loop function can be very easily evaluated from Eq. (fTT^ . for 
instance with Cutkosky rules, giving \mG{{y/s + ie) = — 

In principle Eqs. (USD and m are true only very close to the real axis but, since the 
analytic continuation to general complex plane is unique we can write for a general a/s 


(V^) — Gi {^/s) + i 


A7^^/s' 


(18) 
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with Im(g) > 0. In Eq. ()18|1 one can use for Gj either Eq. (HI or the result of the cutoff 
method. 

One could also have gone to R2 by using Eq. (fT4|l but with the solution of qi with 
Im(gi) < 0, but again one hnds the problem of the multivaluedness of the In functions. 
We have checked, by comparing with the result obtained from Eq. m, that one can use 
Eq. (IT^ as it is written with the prescription of the In explained below Eq. (I14|l and using 
-y/s in the form a + ib, a and b real. 

When looking for poles we will use G{{^) for Re(A/s) < mi + Mi and Gj^{y/s) for 
Re(y^) > mi + Ml. This prescription gives the pole positions and half widths closer to 
those of the corresponding Breit-Wigner forms in the real axis. In this way, when being 
below the lowest threshold, we could also obtain possible pure bound states. 

3.2 SU(3) symmetry breaking scan 

The use of physical vector and pseudoscalar masses, both in the potential Vij as in the 
loop functions, allows for SU(3) symmetry breaking. By following a similar procedure 
to Refs. [17,20], we are going to look for poles starting from the exact SU{3) limit and 
breaking SU{3) symmetry gradually. This can be accomplished by starting using a unique 
vector meson mass, Mq, and pseudoscalar mass, mo, (obtained as an average of the physical 
masses inside each SU{3) multiplet) and changing the masses in the following way: 


M^{x)=M^,+x{M^-M^,) 

mf{x) = ml +x{m'^ - ml), (19) 

with 0 < X < 1. In this way, the masses used are the SU{3) averaged masses for x = 0 
and the physical masses for x = 1. In Fig. HI we show the position of the poles for the 
{S, I) channels (0, 0), (0,1) and (1,1/2), (the only (S', I) channels for which we hnd poles), 
evaluating the loop functions with the subtraction constant method. 

In the exact SU{3) limit (x = 0) we hnd two degenerate octet poles in all the (S',/) 
channels at -y/i = (1168 + iO) MeV, and also a singlet pole for the (0, 0) channel at y/s = 
(1067 + iO) MeV. Hence, our guess at the end of subsection 12.11 that there could be two 
octets and one singlet of dynamically generated resonances gets conhrmed. As we break 
SU{3) gradually, by increasing x, two branches for each (S',/) channels of the octet and 
one for the singlet emerge. The branches end in the physical mass situation, x = 1. Note 
that each branch for the (0, 0) and (0,1) channels has well dehned G—parity while this is 
not the case for (1,1/2) since this corresponds to non-zero strangeness. In the plot we have 
also written for each pole our guess for the correspondence with physical I"*" resonances of 
the PDG. These assignments will be justihed in the detailed discussion of section |3J 

The plot in Fig. HI is instructive but must be interpreted with caution. In principle, the 
real part of the pole position at x = 1 should reflect the mass of the resonance while two 
times the imaginary part should be the width. One may wonder how stable are the results 
with respect to reasonable changes in the regularization scheme. An estimation of the 
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100 



Figure 5: Trajectories of the poles in the scattering amplitudes obtained by increasing the 
SU{3) breaking parameter, x, from 0 to 1. The loop functions are regularized by means of 
the subtraction constant method. 


uncertainties can be done by using the regularization method of fixing G = 0 at the vector 
mass [22], The results are qualitatively similar for the trajectories and the real part for the 
pole positions at a; = 1 differ in the worst of the cases in less than 100 MeV and in most 
cases by less than 30 MeV. The imaginary part of the pole positions also differ in similar 
amounts, but relative to the absolute values of the mass and width of the resonances, the 
differences in the width are more signihcant. Yet, here we must observe that the width 
obtained so far is only a first approximation for the following reasons: 

a) We have not considered other decay channels which are not made of a vector and a 
pseudoscalar. Other channels to which the resonance couples weakly can in practice give a 
sizeable partial decay width because of the large phase space available. For example, this 
would be the case of the 6i(1235) going to four pions, where the phase space is favored 
with respect to our VP channels. While including these extra channels in the coupled 
channels approach goes beyond the scope of the present work, it is important to note that 
the weaker strength of the couplings to these channels makes their repercussion on the real 
part of the pole positions less relevant since there are no restrictions of phase space for the 
real part of the amplitudes. Hence, we might expect small changes in he real part of the 
pole positions from these neglected sources, however, not altering the important fact that 
these poles appear for these quantum numbers. 

b) The second reason is that so far we have not considered the width of the vector 
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mesons in their propagators in the loop functions. For the case of the p and the K* it is 
important to take this into consideration. We shall take this into account in the results 
section and we shall see that this modihes the width of the resonances but only very 
slightly their mass. We will also use another method to account for the width of the vector 
mesons by using the coupling of the resonances to the VP channels and evaluating the 
partial decay widths by means of a convolution with the mass distribution of the vector 
and axial-vector mesons. This will be seen in the next chapter. 


3.3 Couplings and partial decay widths 

The physical interpretation of the poles comes clearer if we realize that close to a pole, 
and if it is not very close to a threshold, the amplitude takes the form of a Breit-Wigner 
structure. By looking at the structure of the covariant amplitude in the Appendix II 
Eq. and the fact that the poles appear in the Vb{l — h) term, close to a pole the 
amplitude has the form (in one channel for simplicity) 


T ~ 


Vh 


fiiy 


pfi pi 




1-bV 

which has the structure of a resonant pole amplitude (see Fig. El) 


( 20 ) 


T ~ 


P2 - Ml 


+ 


piJ.pv 


Ml ) 


( 21 ) 


Generalizing to different channels and considering only the transverse polarizations, the 


Si 


Figure 6: Effective interpretation of the PV PV scattering process close to a pole. 


amplitudes close to a pole in the second Riemann sheet can be expressed as 


p.. ~ 


9i9j 


( 22 ) 


where we have omitted the trivial e ■ e' factor. Hence the factors Pi, which stand as the 
effective coupling of the dynamically generated axial-vector resonances to the channel i, 
can be calculated from the residues of the amplitudes at the complex poles. 

With the values of the couplings we can evaluate the partial decay widths of the axial- 
vector mesons into each VP channels, which reads 


A-^VP — 


\9vi 


SvrMl 


9- 


(23) 
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In the case where there is little phase space for the decay or it takes place due to the 
width of the particles, we fold the expression for the width with the mass distribution of 
the particles as 


1 |•(ArA+2^yl)^ /•(Arv'+2ry)^ 

^A-^VP = ^ / dsA / dsv ^AV\/^) 

J[Ma-2VaY J(Mv-2VvY 


Im 


Im 


_sa- M^ + iMATA) [sv-My + iMvTv 
0 (\/sa - a/sf - Mp)0(yi7 - \/s^)0(a/sf - \/s\. 


(24) 


where 0 is the step function, Tavp = q with q = ^^=X^^‘^{sa, sy, Mp), and Fy 

are the axial and vector mesons total width and are the threshold energies for 

the dominant A and V decay channels. 

In Eq. fl24|l the convolution is done using the masses and total widths of the particles 
taken from the PDG. 


4 Detailed study of the 1+ resonances 

4.1 S = 0,1 = 0 

In Fig. [7| we show the results for the diagonal Ta matrix elements as a function of the 
energy with the loop functions evaluated with the subtraction constant method (we will 
call it a method) (left column), with the prescription of G = 0 at the vector meson mass 
energy (6 method) (central column) and with the cutoff method including the widths of 
the vector mesons in the loops (c method) (right column). This latter method evaluates 
Eq. (IT!?|l performing the integration analytically and the q integration with a cutoff 
Qmax = 1000 MeV substituting Mf —>• Mf — iMiV. 

We do not show the results with the cutoff method without widths in the vector meson 
propagator in the loops since they are very similar to the ones with the a method (actually, 
the subtraction constant has been chosen to agree with the cutoff method). We can clearly 
see the reflexion on the real axis of the poles found in Fig. El in this (0, 0) channel, which 
corresponds to the a method. The a and b methods give qualitative similar results but 
there are signihcant differences at low energies although not so large at higher energies. 
When we include the width in the vector meson propagators of the loops, we obtain a 
broadening of the amplitudes. All these discrepancies give an idea of the uncertainty we 
have in the model. 

In Table El we show the position of the poles obtained with the a method and the 
couplings to the different isospin channels. 

The pole at 919 — il7 MeV, which is found for the channels with negative G—parity, 
may be identihed with the hi(1170) resonance. In the PDG the only decay channel seen so 
far is pn, although we hnd this pole in all the channels allowed by the quantum numbers. 
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Figure 7: Amplitude for S' = 0, / = 0, for the different allowed channels. Left column: 
subtraction constant method in the loop functions; central column: with the prescription 
of G = 0 at the vector meson mass energy; right column: cutoff method considering vector 
meson width in the propagator. 


even if they are kinematically closed. However, in Fig. 0 these other channels are difficult 
to see since the amplitude is dominated in this region by the pvr channel. 

The partial widths obtained for the axial-vector resonance decaying into VP are pre¬ 
sented in Table ini The results that we obtain are compatible with the scarce experimental 
information in the sense that the only decay channel seen is the pvr, which is the only 
relevant in our case. We can see that the widths for this channel calculated with the a and 
b methods are qualitatively similar. It is also instructive to see that with the c method (see 
3rd column of Fig. ED we see the apparent full width of the resonance (around 90 MeV), 
which is in reasonable agreement with the sum of all the partial decay widths in Table El 
for the hi(1170). 

The 1245 — i7 MeV pole can be assigned to the hi (1380) resonance. Note that in the 
PDG the isospin of this resonance is not given, although in the classihcation schemes it is 
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919-il7 

1245 - i7 

1288 - i9 

9i \9i\ 

9i \9i\ 

9i \9i\ 

-V(K*K + K*K) 
(j)T] 

OJT] 

pn 

K*K) 

46 - il3 48 

23 - i28 36 

-3453 + il681 3840 

781 -z498 926 

-3311+ i47 3311 
3020 - i22 3020 

648 - i959 1157 

6147+ M83 6150 

7230 + i9 7230 


Table 5: Pole positions and conplings for S' = 0, / = 0. (All the nnits are MeV.) 



■pexp 

tot 

■pexp 

^ i 

rf(o) 

rf(5) 

hi(1170) ^ K*KPc.c. 

360 ± 40 

— 

0.5 

1.4 

hi(1170) — pn 

// 

seen 

77 

115 

hi(1170) —>• tor] 

// 

— 

0 

0 

hi(1170) ^ (j)r] 

// 

— 

0 

0 

hi (1380) ^ K*K + c.c. 

91 ±30 

seen 

45 

36 

hi(1380) —pvT 

// 

— 

9 

4 

hi (1380) —>• LOT] 

// 

— 

25 

16 

hi (1380) ^ 

// 

— 

0 

0 

/i(1285) ^ k*K + c.c. 

24 ± 1 

not seen 

0.3 

0.2 


Table 6: Partial decay widths for S' = 0, / = 0. (All the nnits are MeV.) 


assnmed to have 1 = 0. Onr assignment, as well as snggested in [22], is also 1 = 0. 

The only experimentally observed decay channel of the hi (1380) is K*K (with G—parity 
negative). In Table El we see onr resnlts for the different partial decay widths. It is difficnit 
to extract concinsions from Table El given the scarce experimental data, bnt we shonld 
stress that the only channel experimentally seen is precisely the dominant one in onr 
calcnlations. In this case the snm of onr partial VP decay channels is compatible with 
the total experimental width. We also see that method c gives a width (aronnd 40 MeV) 
qnalitatively similar to the snm of the partial decay widths in Table El 

The pole at 1288 — iO MeV is below the K*K threshold which is the only allowed 
PV channel for positive G—parity, and therefore the pole appears as a bonnd state of 
K*K + C.C.. Actnally, in Fig. 13 with the a and b methods |Tp goes to inhnity at the pole 
position. It is reasonable to assign this pole to the /i(1285) resonance. The reason why we 
do not get width for this resonance, while in the PDG it is said to be 24 MeV, is becanse 
there are other decay channels different to VP that we are not obvionsly considering, (like 
47 r (33%), rjirn (52%) or KKtt (10%)) . In fact the K*K + c.c. is qnoted in the PDG 
as ’’not seen” which agrees, in a hrst approximation, with the resnlt obtained here as a 
bonnd state. If we evalnate the partial decay width into this channel, with the conpling 
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obtained and with the convolntion of the widths we get abont 0.3 MeV which justihes 
why this decay channel has so far not been seen. This small nnmber is actually rather 
unstable since it can be made even an order of magnitude bigger by enlarging the limits in 
the convolution of Eq. ( 121 . The apparent width from the |Tp plot when considering the 
vector meson width in the propagator of the loop functions is around 20 MeV. With all 
these apparent instabilities in the widths, one should not lose the important point which is 
that, comparatively to the other resonances, the width for this resonance is very small. It 
is worth noting that the /i(1285) couples only to K*K in our theory and with such a large 
strength qualihes strongly as a quasibound K*K or K*K state. The small experimental 
total width indicates a small coupling to other channels which are kinematically open and 
barely change the nature of the resonance as a bound system of K*K. 

When arriving to this point it is important to note that there is experimentally another 
resonance, the /i(1420), which is assigned by the PDG to the nonet. However, in our 
scheme this resonance has no counterpart. This is because, as explained in subsections 12.11 
and l3.2[ the interaction of two octets give one singlet plus two octets (three poles), therefore 
there is no way, with the interaction of one vector and one pseudoscalar mesons, to generate 
dynamically one more pole in this channel and, therefore, there is no room for one more 
resonance like the /i(1420). Therefore, in our scheme, the /i(1420) cannot be considered 
as a dynamically generated resonance from the PV interaction. 

On the other hand, the hi(1170) state that we generate (the original singlet in Fig. El) 
has the same quantum numbers as the isoscalar member of the octet of the 6i(1235), the 
hi(1380), and hence the set of the 6i(1235), hi(1170) and hi(1380), together with the 
strange partners that we will discuss below, agree with the association of a nonet to all 
these states of the PDG. 

4.2 5' = 0, / = 1 

In Fig.IHlwe show the results for the |Tp with the same three methods used in the previous 
subsection. 

We clearly see the reflexions in the real axis of the poles found in the complex plane, 
(shown in Table Ej). 

The pole at 1247 — i28 MeV, which is found for the channels with positive G—parity, 
may be assigned to the 6i(1235) resonance. In the PDG the full width is quoted to be 
142 MeV and the decay channels are quoted as (utt (dominant), prj (seen), dvr (< 50%), 
KKti (< 14%) and 07r (< 1.5%), (See Table |H1). Therefore one could expect that the VP 
decay channels would account for around 50% of the total width. This is compatible with 
what we obtain, as can be seen in Table |H1 Note also that the channel quoted as dominant 
in the PDG is indeed the dominant one in our results. It is also interesting to call the 
attention to the fact that the strongest coupling of the 6i(1235) is to K*K as it was the 
case of the hi (1380), consistently with the fact that they belong to the same SU{3) octet 
in the SU{3) limit. It is interesting to see that this feature remains in spite of the SU{3) 
breaking. The method c provides an apparent total width of about 80 MeV. 

The pole at 1011 — i84 MeV, which is found for the channels with negative G—parity. 
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Figure 8: Amplitude for A = 0, / = 1, for the different channels. Left plot: subtraction 
constant method in the loop functions; central plot: with the prescription of G = 0 at the 
vector meson mass energy; right plot: cutoff method considering vector meson width in 
the propagator. 


\/ Sp 

1011-i84 

1247 - m 

9i \9i\ 

9i \9i\ 

-^{K*K + K*K) 

(j)TT 

UTT 

PI 

pn 

■^{K*K - K*K) 

-3795 + i2330 4453 
1872 - H486 2390 

6172 - i75 6172 

2087 - i385 2122 

-1869 + i300 1893 
-3041 + im 3082 


Table 7: Pole positions and couplings for S' = 0, / = 1. 


could be assigned to the ai(1260) resonance. Note that the values for the mass and width 
quoted for this resonance in the PDG, 1230 ± 40 and 250 — 600 MeV respectively, have a 
large uncertainty. This can justify the discrepancy in the position of the pole we obtain. 
For the partial decays widths there is no good data quoted in the PDG (every one of the 
many decay channels are just quoted as ’’seen”). However, in the detailed explanation in 
the PDG, some experiment gives for pvr 60% and for K*K + c.c. 8 — 15 MeV%, what would 
be in reasonable agreement with what we obtain. The method c provides an apparent 
width of about 300 MeV which is tied to he slow fall down of |Tp at the right of the peak, 
which differ from a typical Breit-Wigner shape. 

It is interesting to note that the octet of the ai(1260) resonance has been considered as 
one of the fundamental helds in chiral theories which deal explicitly with spin-1 helds [25]. 
In this latter work the Lj parameters of [28] have been derived assuming exchange of 
these vector mesons (and scalars). The comment is worth making because the use of the 
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■pexp 
^ tot 

pea;p 

rf(o) 

rf(o) 

6i(1235) ^ K*K + c.c. 

142 ±9 

— 

7 

10 

6i(1235) ^ 07r 

// 

< 1.5% 

12 

13 

6i(1235) ^(UTT 

// 

dominant 

25 

25 

6i(1235) ^ pr] 

// 

seen 

8 

9 

ai(1260) — > pTT 

250 - 600 

seen(60%?) 

106 

156 

ai(1260) — >• K*K + c.c. 

// 

seen{< 10%?) 

6 

11 


Table 8: Partial decay widths for S' = 0, / = 1. (All the units are MeV.) 


chiral interaction between the pseudoscalar and vector mesons used here, together with 
the unitarization procedure, generates dynamically the ai(1260) resonance and, hence, 
introducing it further as an explicit degree of freedom would lead to doublecounting. 

There is another observation worth making which is the fact that there is another 
Oi resonance in the PDG at higher energies, the ai(1640), which we do not generate 
dynamically and is more likely to be a genuine qq meson. At the same time there is 
another observation to make which is the fact that no other bi resonance is quoted in the 
PDG. 

4.3 5 = 1, / = 1/2 

In Fig. El we show the results for the amplitudes with the same three methods considered 
in the other channels. 



Vs (MeV) 


Figure 9: Amplitude for S = 1, I = 1/2^ for the different channels. Left plot: subtraction 
constant method in the loop functions; central plot: with the prescription of G = 0 at the 
vector meson mass energy; right plot: cutoff method considering vector meson width in 
the propagator. 
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In Table E] we show the pole positions and the couplings obtained. The lowest pole 


\/^p 

1112-m 

1216-z4 


9i \9i\ 

9i \9i\ 

(j)K 

1587 - z872 

1811 

1097 - i400 

1168 

ujK 

-1860 + i649 

1970 

-1033 + i375 

1099 

pK 

-1524 + Z1154 

1912 

5274 + i297 

5282 

K*r] 

27 + il55 

157 

3459 - i95 

3460 

K*tt 

4187-i2098 

4683 

340 - i984 

1041 


Table 9: Pole positions and couplings for 5 = 1,1 = 1/2. 


couples strongly to iP*7r, very weakly to K*ri and moderately to the rest of channels. On 
the contrary, the higher pole couples strongly to pK and K*r] while moderately to the 
other channels. 

In Fig. [ini we plot, as an illustrative example, the modulus squared of the scattering 
matrix in the second Riemann sheet for three different channels. The pole structure can 
be clearly seen as well as the relative strength of each channel in the two different poles 
(similar plots can be obtained for the rest of channels). 



Figure 10: |Tp in the second Riemann sheet for (j)K, pK and K*tt channels for 5 = 1, 
I = 1/2. 

In the PDG there are two physical 5 = 1, J = 1/2, resonances with = I"*", which 
are the iFi(1270) and the iFi(1400). Actually these two resonances have been usually 
considered to be a mixture of the Ki members of the I’*"’*' and !■*■“ octets, called Kia 
and Kib respectively. At this point we have difficulty assigning the poles found to these 
resonances. In order to have a feeling of which reasonable assignment to make we study 
the partial decay widths of these two resonances. This is shown in Table In the table 
we have assumed for convolution purposes and phase space the mass of the resonance to be 
1270 MeV and in the last column we show the expected partial decay width coming from 
the poles in Table IHl We can see that the resonance at 1112 MeV couples strongly to the 
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K*Tr and leads to a partial decay width in this channel of 113 MeV (aronnd 150 MeV shonld 
we considered the mass to be 1400 MeV). On the other hand, the 1216 MeV resonance 
conples dominantly to pK. In view of this, one wonld be tempted to assign the 1112 MeV 
pole to the i^i(1400) resonance and the 1216 MeV pole to the i^i(1270) since this is the 
experimental case. Yet, the difference of 300 MeV in the mass in the case of the i^i(1400) 
is not an appealing featnre, since this difference is mnch larger than the i^i(1400) width 
(174 MeV). There is another possible scenario which we hnd more appealing. The detailed 
explanation of the PDG on the determination of the decay widths shows a clear discrepancy 
between different methods of determination of the i7i(1270) width. In particniar, a set 
of experiments nsing K beams leads to mnch larger widths (abont a factor three) than 
another set that nses pion beams. This conld hnd an explanation if one assnmed that the 
experimental i7i(1270) resonance is a snperposition of two resonances that conple with 
different strength to a particniar channel. This is indeed the case for the two poles that we 
hnd, one of them conpling strongly to K*ti and the other one to pK. Diherent experiments 
which favor mechanisms that give a bigger weight to one of these channels would lead to 
very diherent widths of the resonance. This is the case with the recent hndings that the 
A(1405) resonance corresponds to actually two poles and diherent experiments favor one 
or the other pole, leading to diherent visible widths of the resonance [17,18,29]. 

There are some experimental features in the diherent experiments that could be un¬ 
derstood within the two pole structure that we obtain. Indeed, some experiments provide 
a dominance of the i7i(1270) decay into pK [30-32] while other experiments show a clear 
dominance of the K*tt decay mode [33]. In our theoretical approach the dominance of the 
pK decay in a reaction has to be interpreted as a sign that the dynamics of the reaction 
favours the second pole (which couples mostly to pK) and this would have a small width. 
This is indeed the case experimentally, and the total width is about 60 MeV (similar to 
our results). On the other hand, the dominance of the 77 * 7 r decay channel in another 
experiment would have to be interpreted as the dynamics of this reaction favouring the 
coupling to the hrst resonance. In this case the total width would be large and indeed 
this is what is seen in the experiment of [33] where the total width is about 190 MeV (we 
obtain 125 MeV with the a method or 200 MeV with the b method). 

In Ref. [22] a broad bump in the speed plot was associated to the 77i(1400) resonance 
but we hnd no pole in that amplitude in that region. 

It is interesting to see that, when breaking the SU{3) symmetry, the S' = 0 states of the 
two octets in Fig. |Sldo not mix because they have well dehned G—parity (this is opposite 
to the case of the meson baryon interaction [17] where one did not have this constraint). 
However, when we go to S' 7 ^ 0, then the states are no longer eigenstates of G—parity and 
the octets mix. This is a well known case for the Ki axial-vectors where there is much 
discussion about the mixing in the literature^. 

We next proceed to see which is the mixing angle for the two poles in S = 1, / = 1 / 2 , 
that we hnd. We write 

^For a review of the theoretical status on this mixing angle see, for instance, the introduction of Ref. [34]. 
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■pexp 
^ tot 


r‘'*(o) 

rf(6) 

/7i(1270) ^ (j)K 

90 + 20 

— 

(0,0) 

(0.0) 

/7i(1270) ^ a;A: 

// 

11 % 

(5.2) 

(10,0.1) 

/7i(1270) ^ pK 

// 

42% 

(7,52) 

(5,12) 

/7i(1270) ^ K*p 

// 

— 

(0,0) 

(0,0) 

771(1270) ^ K*tt 

// 

16% 

(113,6) 

(194,1) 


Table 10: Partial decay widths for S' = 1, / = 1/2. In the last column, the hrst element is 
the value obtained with the coupling to the lower pole and the second element is the result 
with the coupling to the higher pole. (All the units are MeV.) 


17 ^ 1 ( 1270 ), 1 > = sm9\KiA > + cos9\Kib >, 

1 / 71 ( 1270 ), 2 > = cos9\Kia > —sm9\KiB >, (25) 

where |/7i(1270), 1 > is the state associated to the lower mass pole and |/Pi(1270),2 > to 
the upper mass pole. In Eq. \Kia > corresponds to the SU{3) |8a, 1,1/2, +1/2 > state 
and \Kib > to the |8s, 1,1/2,+1/2 > state, (in the notation \irrep,Y, I, >, with irrep 
being the name of the irreducible representation of SU{3), Y the hypercharge, I the total 
isospin and I 3 the third component of I). (Note that we have chosen = +1/2 although 
the following discussion works obviously for any allowed J3). 

• “ i 

R -. 


Figure 11: Schematic interpretation of the coupling of the dynamically generated axial- 
vector resonance to a PV state. 


The couplings considered so far are (see Fig. HH), by dehnition. 


9. =< fi|(|i(/=l/2,/3 = +l/2) >, 


(26) 


where |/? > is the generic name for the physical resonance associated to either pole. By 
making use of SU{3) Clebsch-Gordan coefficients, we can obtain the relation of the cou¬ 
plings in isospin base to the ones in SU{3) base, giving 
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1112 -ffi4 

1216-74 

927 

9io 

98s 

98a 

—8 + i94 
-346 + ilb 
4259 - i2455 
2317-7928 

316-721 
-406 + 751 
-3783 - 7947 

5208 - 7636 


Table 11: Couplings for S' = 1, / = 1/2, in SU(3) base. 


By using in the left hand side of Eq. the values obtained in Table E we get the results 
shown in Table im 

As mentioned after Eq. (unD, the singlet component of the </> and uj (see Eq. dH}) does 
not lead to an interaction term, and this is already assumed in Eq. (EH) where only the 
octet component of the (j) is taken. We could as well have taken instead of in Eq. 
the coupling —\/2g^K and in the exact SU{3) limit this would lead to the same results. 
When SU (3) is broken, the use of either coupling in the Eqs. (I77|) will give us an idea 
about the uncertainties in this SU{3) decomposition. 

In Table [TT] we can see that the couplings to the 27—plet and the decuplet are almost 
negligible in comparison to the two octets, in fact, they are compatible with zero within 
the uncertainties we can assume in our model. This means that the two poles found are 
essentially a mixing of the antisymmetric and symmetric octets. The above reasoning is 
more than qualitative since it allows us to quantify the weight of the SU (3) components of 
the resonance and therefore the mixing angle dehned in Eq. (ESI). In order to obtain this 
mixing angle let us write Eq. (ESI) in a more generic way: 

\R >= a27|27, R > +aio|10, R > Tass |8s5 R > +®8a R > (28) 

where we have used the short notation |27, i? > = 127,1,1/2,+1/2,7? > and so on for 
the rest. The aj coefficients must satisfy 1%^ = 1- To obtain the mixing angle 6 we 
have to evaluate as^ and as^ since, by dehnition, they are as^ = cos 9 and as^ = sinT' for 
|A:i(1270),1 >. 

On the other hand, up to lowest order in SU{3) breaking, we have 


gss = < ^l^|8s .RN8^ >= a), 

98a = < as, < 8a, R\t\Sa >= aljsa b). (29) 

Therefore, had we know fs,, ^8,, we would know as,, as^, by using gg^, gsay from Table ITTl 
To evaluate fs, and fs, we can go to the SU{3) limit, which by virtue of Eq. (fTT|l should 
be the same. In this limit the two octet poles are in the same place. The couplings ts,, 
tga, are easily evaluated since the matrix elements of the potential in the SU{3) base are 
already known in Eq. EH). Hence, by performing the Bethe-Salpeter resummation in this 
base, we readily obtain the poles corresponding to each SU{3) irreducible representations 
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and the couplings. The amplitudes behave in the pole as 


< 8.|r|8, >= 



( 30 ) 


and analogously for 8a and singlet, the only channels where there are poles. We hnd 

= ^8a = 5568 MeV. 

In the case of physical masses, a; = 1, using the same values for and fg^ that we have 
obtained in the SU{3) limit and the couplings gs^, of Table ITTl we obtain the following 
results: for the pole 1112 — i64 MeV we get from Eq. (j21r) e ~ 28°, and from Eq. (|2nb) 
0 ~ 27°. Should we have used —\/2g^K instead of g^x in the calculations we would get 
0 ~ 22° from Eq. (I29h l and 0 ~ 34° from Eq. (&)• These discrepancies give an idea of 
the uncertainties in our calculations. 

For the pole 1216 — i4 MeV, we get from Eq. ipilk l 44°, and from Eq. (jSSb) 20°. And 
with g^x 43° and 11° respectively, which implies that we have a larger uncertainty in this 
case. In fact, for the pole 1112 — i64 MeV we get X) kiP = 0-98 by using g^x and 1.17 by 
using g^jx^ and for the pole 1216 — iA MeV, 1.38 and 1.43. The deviation of these numbers 
with respect to 1 gives an idea of the uncertainties. 

I summary, we get a mixing angle of the order of 30° with an uncertainty of about 40%. 
We should however refrain from making any association of the mixing angle found here 
with the mixing angle of Kia, Kib, mentioned in the literature since this latter one is used 
to produce the 77i(1270) and 771(1400) resonances and here we are mixing Kia, Kib, to 
produce |i7i(1270),l > and |i7i(1270),2 >. 


Along this work we have devoted some attention to uncertainties in the theory. We 
would like to call the attention to another source of uncertainties which would affect all 
channels. So far we have used only the pion decay constant, f = f.^ = 92 MeV, in 
the Lagrangian. We could have used different values for /^, fx, fri- In order to estimate 
uncertainties from this source we have proceeded as in [13] and taken an average / = l.lb/jr. 
We hnd that all the poles obtained so far still appear but the pole positions are somewhat 
changed. The trajectories of the poles in the complex plane are essentially the same 
but with real parts shifted about 50 MeV to higher energies (using the same subtraction 
constant). The imaginary parts are accordingly increased since there is now more phase 
space. The most signihcant changes in the imaginary part of the pole positions are for the 
poles associated to the hi(1170) (100% increase), ai(1260) (40% increase) and i7i(1270; 1) 
(50% increase) and small changes in the rest. These increases are mostly tied to the 
increased phase space. However, we have so far evaluated the partial decay widths of the 
resonances, using the couplings obtained and the physical masses of the resonances. We 
have checked that the couplings of the resonances barely change with the use of the new / 
constant (less than 10% change in general and less than 5% in the dominant channels) and 
hence, within the uncertainties discussed along the work, the results and the conclusions 
are unchanged from this new source of uncertainty. 

The overall conclusion of all these tests is that the existence of the poles and their 
basic properties are very solid and not contingent to the difference sources of uncertainties 
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discussed along the paper. 


5 Conclusions 

In this work we have done a systematic search for possible I"*"”, dynamically 

generated resonances through the interaction of vector and pseudoscalar mesons. The 
starting point has been a chiral Lagrangian which, by expanding up to two pseudoscalar 
helds, leads to a Weinberg-Tomozawa term which accounts for the two-vectors and two- 
pseudoscalar meson interaction. From this Lagrangian we have argued, by going to the 
SU{3) limit that from the interaction of the two octets we could expect attraction for one 
singlet and two octets. After that, we have implemented unitarity in coupled channels to 
account for the resummation of loops in L = 0. This resummation has been accounted 
for by means of a Bethe-Salpeter like equation in coupled channels, where some subtleties 
in the evaluation of the loop functions coming from the use of vector mesons have been 
discussed. The regularization of the loops has been done with dimensional regularization 
by means of a subtraction constant hxed to agree with the numerical result obtained 
performing the integration with a cutoff of natural size. We also compare our method with 
another one where the loop functions are hxed to zero when ^/s is equal to the vector 
meson mass. This served to have an idea of uncertainties in the theory. 

We have looked for poles of the scattering amplitudes in the second Riemann sheet. In 
the SU{3) symmetric case, considered by taking equal masses for the vectors and equal 
masses for the pseudoscalars, we found two poles in the same position corresponding to 
two degenerate octets and one pole corresponding to the singlet. As SU{3) symmetry is 
gradually broken, the two degenerate poles split apart in trajectories (different for each 
{S, I) channel) ending in the physical situation, when the physical masses are used. 

By evaluating the residues of the amplitudes in the poles, we have obtained the cou¬ 
plings of the axial-vector dynamically generated resonances to each VP channel, which has 
allowed us to evaluate the different partial decay widths. In view of the information sup¬ 
plied by the pole positions, couplings and partial widths, we have done a correspondence of 
the poles found to the hi(1170), /ii(1380), /i(1285), &i(1235), ai(1260) and A'i(1270). For 
this latter case we found actually two poles coupled strongly to K*tt and Kp respectively 
which also differed appreciably in the width. We suggested that different experiments give 
different weight to each of these resonances and this could explain the discrepancies in 
the widths obtained in different reactions. This would also explain the correlation found 
between experiments finding a dominance of the pK decay (which produce a small width) 
and those finding a dominance of the K*ti decay mode (which produce a large width). We 
also evaluated the couplings of our two iFi(1270) states, hnding a reasonable mixing of 
around 30°. 

The only axial-vector resonances for which we do not hnd poles are the /i(1420) and the 
A'i(1400) for different reasons: from the interaction of two octets one can only generate 
one singlet and two octets, therefore there is no room for more poles apart from those 
discussed above. Actually we only found a pole with the /i quantum numbers which suits 
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better to the /i(1285). In the two poles that we hnd in the S = 1, I = 1/2^ channel there 
are no clues to identify one of them to the it'i(1400) resonance. 

The conclusions reached in this paper about the dynamical nature of these axial-vector 
mesons should have experimental repercussions in the sense that, with the information 
obtained in the present work, one can make predictions for production of these resonances 
in different reactions, which are amenable of experimental search. This has been the case 
for other dynamically generated resonances and we hope the present paper encourages 
work in this direction. 
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Appendix I: S-wave and on-shell VP ^ V'P' tree level 
amplitude 

Let us evaluate the s-wave projection of the tree level amplitude for the process V{q)P{p) 
V{q')P{p'). 

After performing the SU{3) trace in Eq. [7|we get the general expression 

c = -^Cij{d'^v^d,pv'^p' - d'^v^pv'^d.p' - v^d.Pd'^v'^^p' + v^Pd'^v'^^d.p'), (3i) 

where now lA, P and P' are the meson helds, (not the SU{3) matrices of mesons). 
Eq. (pn|l leads to the following amplitude: 

V = -^Vj{p + p'){q + P = - u), (32) 

where s = {p + q)^ = {p' -|- q'Y and u = {p' — qY = {q' — pY are the usual Mandelstam 
kinematical variables. 

The partial wave expansion of the amplitude can be written as 

T = J2{2l + l)fi{s)Pi{x), (33) 

with X = cos6', 6 the center of mass scattering angle and Pi are the Legendre polynomials. 
Hence, the s-wave projection of the scattering amplitude is 

fi=o{s) = ^ ^ T{s, t{x'),u{x'))Pi=o{x')dx' (34) 


27 





The I = 0 partial wave is what we will call the potential Vij. 

Expressing u in terms of x and taking the momenta on-shell, we have u = -|- 

— 2E{p')E{q) — 2|i?’| l^x, and the term proportional to x vanishes when performing the 
integration, and thus the s-wave projection of the amplitude reads 

VtlM = -JwC.iU - + 2E(p')E{g)), (35) 

with 

E(p') = + ; E{q) = ^{s + - m^). (36) 

2y/ s 2 Y s 

Altogether, the tree level on-shell and s-wave amplitude is 






r/2 


3s - (M^ +m^ + M'^ + m'A - - m^)(M- - m 

s 


/ 2 \ 


(37) 


Appendix II: Vector mesons polarization in the re¬ 
summation of the loops 

In this Appendix we are going to explain some subtleties about the treatment of the 
polarization vectors of the vector mesons in the resummation of the loops in the Bethe- 
Salpeter equation. 

First let us show that, in order to look for poles, we can deal only with the transverse 
vector-meson polarization modes. Let us consider the case with only one channel (the 
generalization to coupled channels is straightforward). Let us consider the loop diagram 

of Fig. 1121 


P-q 



q 

Figure 12: Vector-Pseudoscalar loop. 


The Feynman rule is 


t = Ve^i 


d‘^q 1 1 

(27r)"^ (P — qy — w? + ie q^ — -|- ie 




+ 


Ip) 


Vey 


(38) 


where we have taken for the tree level t = Ve ■ e'. 

Note that we have factorized on-shell the V functions which, as shown in [8] using 
the Bethe-Salpeter equation or in [10] using the N/D method, leads to a well dehned 
renormalization scheme. 
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Since the integral only depends on an external momentum P, the result of the integral 
must be of the form 

+ c'P^P^ (39) 

which can be written as 

+ cC^'^) (40) 

with and where we have factorized one of the vertices, V. 

The B and C tensors are idempotent and orthogonal in the sense that 

= P'^^, 

(41) 


Therefore in the iteration of the loops we have two independent series that sum up to 


Vb ( p/^P^\ Vc P^P^ 

1-6 V ^1-c F 2 


(42) 


In the center of mass frame we have P^ = (-^/s, 0, 0, 0). The transverse vector mesons 
(which have only space components) only contribute to the left term of Eq. ()42j) . while the 
longitudinal vector mesons (which has also time component in the polarization vectors) 
contribute to both terms of Eq. (m. 

Therefore, when looking for poles of the T matrix, since if there are poles they should 
exist for all modes and the second term of Eq. does not contribute for transverse 
modes, the poles can only be in Vb/ {1 — b) which is common for all modes. In practice we 
hnd that c is one order of magnitude smaller than b and of opposite sign and the second 
term of Eq. dH leads indeed to no poles. 

Therefore, in order to look for poles we can only consider transverse polarization vectors. 
In this case the tree level VP potential can be written a.st = Ve-e' = -Ve-P, and Eq. 
reads 


t = —Vei i 


_A_1_ 1 (-V) ‘ 

(27r)'^ (P — qY — Tn? + ie + ie \ / 


(43) 


The term qiqj can be replaced by l/3q‘^6ij. In addition, as shown in Ref. [35], 
term can be taken on-shell {Yon — l/(4s)A(s, M^, m^)) since the off-shell part q^ 
can be reabsorbed into the renormalization of the couplings. Therefore Eq. (IlHl) 
expressed as. 


t = V‘^eieiG{l -|- 


^ Y'on 
3 M2 


5 


the g ^ 

- T2 
y ON 

can be 
(44) 


with 


G'(^/i) = i 


Pq 1 1 

(27r)^ (P — g)2 — m2 -|- ie g2 _ ]\g'2 _|_ 


(45) 
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Therefore the series of Fig. ^ would give 


T ——Veiei + (—+ 


1 Qon 
3 M2 


) + {—V)^eieiG^{l + 


Iq 


ON\2 


3 M2 


r + 


(46) 


which sums up to 


T = 


M 


1 + VGil + l^-^) 


e ■ e 


(47) 


Eq. (|T7jl can be easily generalized to more than one channel giving a Bethe-Salpeter 
like coupled channel equation 


T = [1 + VG]-\-V)€-€', (48) 

where G = (7(1 + 3 ]^) is a diagonal matrix with the /—th element, Gi, being the two 

‘ 1 5^ 

meson loop function containing a vector and a pseudoscalar meson. The term is small 

and has no much repercussion in the final results. By comparing Eq. ()48|1 and Eq. (11^ we 
can see that b = —VG. 
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